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ABSTRACT 

Aims. We study the dwarf spheroidal galaxies in the nearby M81 group in order to construct their photometric metallicity distribu- 
tions and to investigate the potential presence of population gradients. We select all the dwarf spheroidals with available Hubble Space 
Telescope / Advanced Camera for Surveys archival observations, nine in total. 

Methods. We interpolate isochrones so as to assign a photometric metallicity to each star within a selection box in the color-magnitude 
diagram of each dwarf galaxy. We assume that the dwarf spheroidals contain mainly an old stellar population. In order to search for 
metallicity gradients, we examine the spatial distribution of two stellar populations that we separate according to their metallicities. 
Results. As a result, we present the photometric metallicity distribution functions, the cumulative histograms and smoothed density 
maps of the metal-poor and metal-rich stars as well as of the intermediate-age stars. 

Conclusions. From our photometric data we find that all the dwarf spheroidals show a wide range in metallicities, with mean val- 
ues that are typical for old and metal-poor systems, with the exception of one dwarf spheroidal, namely IKN. Some of our dwarf 
spheroidals exhibit characteristics of transition-type dwarfs. Compared to the Local Group transition type dwarfs, the M81 group 
ones appear to have mean metallicity values slightly more metal-rich at a given luminosity. All the dwarf spheroidals considered here 
appear to exhibit either population gradients or spatial variations in the centroids of their metal-poor and metal-rich population. In 
addition, there are luminous AGB stars detected in all of them with spatial distributions suggesting that they are well mixed with the 
old stars. 

Key words. Galaxies: dwarf - Galaxies: groups: individual: M 81 group - Galaxies: stellar content - Galaxies: structure - Galaxies: 
interactions 



1. Introduction 

The dwarf galaxies within our Local Group (LG) have been the 
subject of intensive spectroscopic and photometric observations 
in different wavelength regimes and thus are well studied ob- 
jects. Their study has been facilitated by the proximity of these 
I dwarfs so that individual stars may even be resolved down to the 
, main sequence turn-off, depending on their distance. Thus, ex- 
tending the studies to dwarf galaxies in nearby groups with dif- 
ferent environment and comparing their properties are of great 
importance in order to understand the main drivers of their evo- 
' lution. In addition, the derived properties can provide a way to 
constrain models of galaxy formation as well as chemical evolu- 
tionary models 

In this respect, the M 8 1 group is an interesting target: despite 
several differences, it bears close resemblance to our LG. The 
similarity of the M 8 1 group to our LG lies in its binary structure 
(Karachentsev et al. 2002 ), while its difference is mainly due to 
the recent interactions between its dominant consituents as re- 
vealed by the formation of tidal tails and bridges detected in HI 
observations (Appleton, Davies & Stephenson il981: Yun, Ho & 
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Lo'1994i With a mean distance of -3.7 Mpc (Karachentsev et 
al. 2002 ), the M81 group is one of the nearest groups to our 
own LG. It consists of about 40 dwarfs of both early-type and 
late-type, with the addition of 12 recently discovered dwarf can- 
didates (Chiboucas, Karachentsev & Tullv l2009] l. 

Here we focus on the dwarf spheroidal galaxies (dSphs) 
in the M81 group with available Hubble Space Telescope 
(HST) / Advanced Camera for Surveys (ACS) archival data. The 
dSphs are objects with low surface brightness and poor in gas 
content. For a summary of their properties we refer to Grebel, 
Gallagher & Harbeck (2003; and references therein). We use 
their color-magnitude diagrams (CMDs) to derive the photo- 
metric metallicity distribution functions (MDFs) and search for 
the potential presence of population gradients in the M 8 1 group 
dSphs. 

The use of the CMD to infer the star formation histories and 
MDFs is a very powerful tool. For nearby groups at distances, 
where individual red giants are beyond the reach of spectroscopy 
even with 8-10 m class telescopes, CMDs are the best means 
to derive evolutionary histories. With the use of HST observa- 
tions of adequate depth, the upper part of the red giant branch 
(RGB) can be resolved into single stars. Many studies have de- 
rived the photometric MDFs of distant LG dSphs (for example 
Cetus by Saraje dini et al. l2002l And VI and And VII by Grebel 
& Guhathakurta 1 1 9991 ) from their CMDs. A similar work to de- 
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rive the photometric MDFs for dwarf galaxies in nearby groups 
has not been done so far. 

The search for radial population gradients in LG dwarf 
galaxies has been favoured by the fact that the resolved stellar 
populations reach the horizontal branch or extend even below the 
main-sequence turn-off depending on the distance of the dwarf, 
permitting one to use a variety of different stellar tracers. There 
are several studies for population gradients in the LG dwarfs and 
as an example of such studies we refer to the work done by 
Hurley-Keller, Mateo & Grebel ( fT999l l. Harbeck et al. ("SOOB, 
BattagHa et al. (120061) (photometric) and Tolstoy et al. ( 200^1, 
Koch et al. ( 2006) (spectroscopic). There is not any study so far 
searching for population gradients in nearby group dwarf galax- 
ies. 

This paper is structured as follows. In §2 we present the ob- 
servations, in §3 we show our results, in §4 we discuss our main 
findings and in §5 we present our conclusions. 

2. Observations 

We use HST/ACS archival data that were retrieved through 
the Multimission Archive at STScI (MAST). The details of the 
datasets used are listed in Table[T] where the columns show: (1) 
the galaxy name, (2) and (3) equatorial coordinates of the field 
centers (J2000.0), (4) the number of the Program ID and the PI, 
(5) the ACS/WFC filters used, and (6) the total exposure time 
for each filter. 

The data reduction was carried out using Dolphot, a modi- 
fied version of the HSTphot photometry package (Dolphin 2000 ) 
developed specifically for ACS point source photometry. The re- 
duction steps followed are the ones described in the ACS module 
of the Dolphot manual. In the Dolphot output photometric cata- 
logue, only objects with S /N > 5 and "type" equal to 1, which 
means "good stars", were allowed to enter the final photomet- 
ric catalogue. The "type" is a Dolphot parameter that is used to 
distinguish objects that are classified as "good stars", "elongated 
objects", "too sharp objects" and so on. After this first selection, 
quality cuts were applied so as to further clean the photomet- 
ric catalogue. These cuts were based on the distributions of the 
sharpness and crowding parameters, as suggested in the Dolphot 
manual and also in Williams et al. (2009 ). Guided by these dis- 
tributions, we use for the sharpness parameter the restriction of 
\sliarpness filler + sharpnessff;i4w\ < a, where (1.0 < c < 1.5) 
depending on the dSph, and for the crowding parameter the re- 
qniiement (Crowding filler + Crow(i/ng/r8i4vy) < 1.0, where "fil- 
ter" corresponds to either the F606W or the F475W filter. These 
selections were made so as to ensure that only stellar objects 
have entered our final photometric catalogue. The number of 
stars recovered after applying all the photometric selection cri- 
teria are listed in Table|2] column (3). 

The photometry obtained with Dolphot provides magnitudes 
in both the ACS / WFC and the Landolt UB VRI photometric sys- 
tems using the transformations provided by Sirianni et al. (2005 ) 
for the UB VRI system. In the analysis presented throughout this 
work, we chose to use the ACS / WFC filter system. Therefore, 
if we use data from the literature computed in the UBVRI pho- 
tometric system, we transform them to the ACS /WFC system. 
There are two cases where this is necessary. The first case is 
the extinction. The galactic foreground extinction in the V-band 
and I-band, Aj and Ay, taken from Schlegel, Finkbeiner & Davis 
([1998) through NED, are transformed into the ACS/WFC sys- 
tem. For the transformation, we use the coiTesponding extinction 
ratios A{P) / E{B - V) for a G2 star, where A{P) corresponds to 
the extinctions in the filters F814W and F606W (or F475W), 



which are provided by Sirianni et al. (120051 their Table 14). We 
note that the assumption of a largely color-independent redden- 
ing for the RGB is justified since theoretical models indicate 
that the expected eff'ect of the change of color accross the RGB 
amounts to at most 0.01 in E{V - I) for our data (see Grebel 
& Roberts ri995l l. We multiply these extinction ratios with the 
E(B - V), in order to finally get the extinctions in the ACS fil- 
ters. The transformed values, AfSMw and Af(,o6w (or Af4T5w), 
are listed in Table |2] columns (6) and (7) respectively. 

The second case is the I-band tip of the RGB (TRGB), 
which we transform to the F814W-band TRGB in the following 
way. As already mentioned, Dolphot provides the magnitudes 
both in the instrumental ACS/WFC system and in the trans- 
formed UBVRI. Thus, in the range of magnitudes near the I- 
band TRGB, we compute the difference in magnitudes between 
the I-band and F814W-band. This difference is 0.01 mag for all 
the dSphs except for DD044 and DDO 78, where the difference 
is -0.015 mag. The F814W-band TRGB is then equal to the sum 
of this difference and the I-band TRGB. We confirm further more 
this approach of estimating the F814W-band TRGB by applying 
a Sobel-filtering technique to the luminosity function of some 
of the dSphs (Lee, Freedman & Madore 1993; Sakai, Madore 
& Freedman , 1996) estimating the location of the F814W-band 
TRGB. We find that these approaches give values that are in 
good agreement, with a mean difference between them of the 
order of 0.05 mag. 

We list the global properties of the present dSph sample in 
Table |2] The columns contain the following information: (1) 
the galaxy name, (2) the galaxy type, (3) the number of stars 
detected after applying all the photometric selection criteria, 
(4) the visual absolute magnitude My of each galaxy adopted 
from Karachentsev et al. (1200011200111, A lonso-Garcia, Mateo & 
Aparicio (2006), Georgiev et al. ("2005^, (5) the I-band TRGB 
adopted from Karachentsev et al. (2006 200T] l2000l [T999] l. 
(6) and (7) the foreground extinction derived by us for the 
ACS/WFC filters F814W, F606W and F475W, as described in 
Sec. |2] (8) the true distance moduli adopted from Karachentsev 
et al. (l2U56ll200n 120001 [T999] l. (9) the deprojected distance of 
the dSphs from the M 8 1 galaxy, R, adopted from Karachentsev 
et al. (2002 ), (10) the effective radius, rg// adopted from Sharina 
et al. (2008.) and Karachentseva et al. (1987). The dSphs in 
Table|2]are sorted according to their My value. 

Finally, the pixel scale of the ACS/WFC is 0.05" with a 
field of view of 202" x 202"or 4 096 x 4 096 pixels. Thus for the 
mean distance of ~3.7 Mpc of the M 81 group (Karachentsev et 
al. '2002) this field of view coiTesponds to 3.6 kpc x 3.6 kpc, or 
simply 1 pixel corresponds to roughly 1 pc. 

3. Results 

3. 1 . Color-Magnitude Diagrams 

We show the CMDs of the nine dSphs in Fig.[T] where we note 
the difference in the x-axis. The proximity of the M81 group 
and the depth of the observations allow us to resolve the upper 
part of the RGB into individual stars. The most prominent fea- 
ture seen in our CMDs is the RGB. We note the presence of stars 
above the TRGB, which is indicated in Fig.[T]with a dashed line. 
These stars are most likely luminous AGB stars, which indicate 
the presence of stellar populations in an age range from 1 Gyr up 
to less than 10 Gyr. In addition, some of the dSphs appear to have 
bluer stars that probably belong to a younger main sequence. The 
dwarfs in our sample are classified as dSphs (Karachentsev et 
al. l2004l l. with the exception of KDG 6 1 , KDG 64, DDO 7 1 , and 
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Table 1. Log of Observations. 



Galaxy 


RA(J2 000.0) 


Dec (J2 000.0) 


Program ID / PI 


ACS/WFC Filters 


Exposure time 
(s) 
(6) 


(1) 


(2) 


(3) 


(4) 


(5) 


KDG 61 
KDG 64 
DD0 71 
F12D1 
F6D1 


09 57 03.10 

10 07 01.90 
10 05 06.40 
09 50 10.50 
09 45 10.00 


+68 35 31.0 
+67 49 39.0 
+66 33 32.0 
+67 30 24.0 
+68 45 54.0 


GO 9 884 / Armandroff 


F606W/F814W 


8 600/9000 


HS 117 
IKN 


10 21 25.20 
10 08 05.90 


+71 06 51.0 
+68 23 57.0 


SNAP 9 771 / Karachentsev 


F606W/F814W 


1200/900 


DD078 
DD044 


10 26 28.00 
07 34 11.50 


+67 39 35.0 
+66 52 47.0 


GO 10 915/Dalcanton 


F475W/F814W 


2 274/2292 
2 361/2430 



Note.- Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes and arcseconds. 
Table 2. Global Properties (see text for references). 



Galaxy 


Type 


A'. 


Mv 


Itrgb 




A a 
^F606W^ 


(m - M)o 


R 


reff 








(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(kpc) 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


KDG61 


dirr / dSph 


53 543 


-13.87 


23.86 + 0.15 


0.131 


0.202 


27.78 + 0.15 


44 


48 


KDG 64 


dirr / dSph 


38 012 


-13.43 


23.90 + 0.15 


0.099 


0.152 


27.84 + 0.15 


126 


28 


DD071 


dIrr / dSph 


37 291 


-13.22 


23.83 + 0.15 


0.173 


0.267 


27.72 + 0.15 


211 


59 


F12D1 


dSph 


39 519 


-12.84 


23.95 + 0.15 


0.263 


0.404 


27.71+0.15 


181 


31 


DD078 


dSph 


21073 


-12.83 


23.85 + 0.15 


0.040 


0.079 


27.85 +0.15 


223 


38 


DD044 


dSph 


19 357 


-12.56 


23.55 + 0.15 


0.075 


0.149 


27.52 + 0.15 


901 


28 


IKN 


dSph 


14 600 


-11.51 


23.94 + 0.15 


0.111 


0.171 


27.87 + 0.18 


110 




F6D1 


dSph 


14 260 


-11.46 


23.77 + 0.14 


0.144 


0.222 


27.66 + 0.17 


218 


32 


HS 117 


dirr / dSph 


4 596 


-11.31 


24.16 + 0.15 


0.210 


0.323 


27.99 + 0.18 


204 


29 



" or Af„5w in the case of DDO 44 and DDO 78 



HS 117 which are classified as transition-types (dIrr/dSph) as 
they have detectable HI content (Huchtmeier & Skillman ll998l 
Boyce et al. .2001) or Ha emission (Karachentsev & Kaisin 
I20U71 Karachentsev et al. l2006] l. 

3.2. Photometric Metallicity Distribution Functions 

We show the photometric MDFs for the nine dSphs in 
Fig.|2] These are constructed using linear interpolation between 
Dartmouth isochrones (Dotter et al. 2008 ) with a fixed age of 
12.5 Gyr We use Dartmouth isochrones, since they give the best 
simultaneous fit to the full stellar distribution within a CMD as 
demonstrated by e.g. Glatt et al. (I2008al I2008bl l. We chose the 
fixed age of 12.5 Gyr since the RGB in these dSphs may be as- 
sumed to consist of mainly old stars in an age range of about 
10 Gyr to 13 Gyr. The assumption of an old isochrone is also 
justified by the omnipresence of old stellar populations in all of 
the LG dSphs (Grebel|200lJ Grebel & Gallagher |200l and by 
the comparatively small number of luminous AGB stars above 
the TRGB. 

The choice of 12.5 Gyr is an assumption we are making in 
order to estimate the MDFs for each dwarf, while the choice of 
another age in the above range would not considerably affect our 
results. Indeed, the colors of the stars on the RGB are mostly af- 
fected by metallicity differences rather than age differences (see 
for example Caldwell et al. IT9981 Harr is, Haiiis & Poole [19991 
(their Fig. 6); Frayn & Gilmore "SOOT (their Fig. 2)). Thus the 
observed spread in the RGB color is likely caused by a metallic- 
ity spread rather than an age spread, justifying our choice of an 
constant age isochrone. 



The isochrone metallicities we used range from -2.50 dex 
to -0.50 dex with a step of 0.05 dex. The isochrone step we 
used is chosen such that it is smaller than the photometric errors. 
Representative photometric error bars are indicated with crosses 
in Fig. [T] To account for the influence of crowding and the pho- 
tometric quality, we conducted artificial star tests. For that pur- 
pose, we used the utilities provided and as described in Dolphot. 

Given the lack of any spectroscopic information, this method 
of deriving the MDFs is in general fairly accurate as discussed 
in Frayn & Gilmore (I2002I I. In practice, we interpolate between 
the two closest isochrones bracketing the color of a star, in order 
to find the metallicity of that star 

We only select stars within a box plausibly containing stars 
on the upper RGB to construct the galaxies' MDF The bright 
magnitude limit of the box is chosen to exclude the stars brighter 
than the TRGB which belong mainly to the luminous AGB 
phase. The faint magnitude limit of the box is chosen to fulfil 
the requirement that the formal error in the derived [Fe/H] is less 
than 0.15 dex, or 0.2 dex in the case of IKN and HS 117 when 
the photometric errors are taken into account. We employ a dif- 
ferent selection criterion in the case of IKN and HS 1 17 in order 
to have a significant number of stars in their sample as compared 
to the one of the remaining dSphs. The selection criterion based 
on the metallicity formal error depends on the depth of the ob- 
servations. In our data sample we distinguish three categories, 
from now on referred to as depth categories. The first depth cat- 
egory contains KDG61, KDG 64, DD071, F6D1 and F12D1. 
The second depth category contains DDO 44 and DDO 78. The 
third depth category contains IKN and HS 117. Each depth cat- 
egory contains those dSphs that belong to the same Program ID 
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Fig. 1. Color-magnitude diagrams for the nine dSphs. The horizontal dashed lines show the location of the TRGB. The crosses on 
the left hand side correspond to the photometric eiTors as derived from artificial star tests. The boxes enclose the stars for which we 
derive the photometric metallicities. 



and thus have the same filters and roughly the same exposure 
times, as listed in Table[T] 

In order to estimate the faint magnitude limit for each dSph's 
RGB box as a function of the eiTor in [Fe/H], we proceed as fol- 
lows. We extend the faint limit of the bounding box to a mag- 
nitude limit of 26 in F814W for all the dSphs. We compute the 
[Fe/H] for all the stars within each dSph's box, as well as the 
corresponding errors in metallicity. We show the derived mean 
errors in metallicity versus the F814W-band magnitude in Fig. [3] 
for KDG61, DD044, IKN and HS 117, which are chosen here 
as representative examples of the three depth categories. In the 
case of IKN and HS 117, which belong to the third depth cat- 
egory, we show the corresponding plots for both since the re- 



quirement of 0.20 dex leads to slightly different faint limits of 
the RGB box. Based on these plots and on the metallicity re- 
quirements, we choose 25 and 24.5 mag as faint limit for the first 
and second depth category, while in the case of IKN and HS 117 
we choose 24.4 and 24.6 mag, respectively. The choice of these 
limits corresponds to an error in color of less than 0.02 mag for 
the first depth category and less than 0.07 mag for the remaining 
two depth categories. We note that the difference in the eiTor in 
color is due to the different exposure times for each dSph data 
set, which are listed in the column (6) of Table[T] The RGB boxes 
used in each dSph are drawn in Fig.[T] 

In Fig. m left panel, we plot the RGB box used in the case 
of the dwarf KDG 64 as well as a subset of the isochrones used 
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Fig. 2. Metallicity distribution functions for the nine dSphs sorted by their absolute My magnitude, from top to bottom and from left 
to right. The solid lines show the metallicity distribution convolved with the errors in metallicity. The dashed lines show the fitted 
gaussian distributions. The error bars in the upper right corner, or upper left in the case of DD044, show the 1 cr spread for the 
weighted mean metallicity we derive from our data. Note the different scaling of the individual y-axes. 
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Fig. 3. Mean error in metallicity versus the F8 14W-band magni- 
tude. The circles indicate the mean error in metallicity in each 
magnitude bin, while the error bars indicate the standard devia- 
tion of the errors. 



for the interpolation method, here ranging from -2.50 dex to 
-0.80 dex. The step remains 0.05 dex. The grid of the theoretical 
isochrones we use is fine enough that the spacing between them 
is kept nearly vertical. We coiTect the magnitudes and colors of 
the theoretical isochrones for foreground Galactic extinction and 
for the distance modulus of each dSph. The A^guiv and A/7606w 
(or Af475vi' in the case of DD044 and DD078) that we calculate 
and the true distance moduli are listed in Table |2l columns (6) 
and (7) for the extinction and (8) for the distance moduli. The I- 
band TRGB values shown in column (5) of the same table were 
used to compute the F814W-band TRGB values, as explained 
already in Sec.|2] 

We note in Fig. |4] the presence of stars within the RGB box 
bluewards of the most metal-poor isochrone available from the 
Dartmouth set of isochrones. These stars are not used to con- 
struct the MDF. The existence of such stars may indicate the 
presence of more metal-poor RGB stars or of old AGB stars, 
with ages typically greater or equal to 10 Gyr. Such old AGB 
stars that have the same luminosity as RGB stars were also 
noted, for example, by Harris, HaiTis & Poole ( 1999) while con- 
structing the MDF for stars in a halo field of the giant ellipti- 
cal NGC5128. It is expected that at most 22 % of the stars in 
the RGB selection box, and within 1 mag below the TRGB, 
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F606W-F814W [Fe/H] 



Fig. 4. Left panel: Color-magnitude diagram for KDG64 
zoomed in the RGB part to show the stars selected in the 
box, shown with the dashed line, for which we compute their 
photometric metallicities. In the same figure we overplot with 
solid lines a subset of the isochrones used for the interpolation 
method, with metallicities ranging from -2.50 dex to -0.80 dex. 
Right panel: The error in metallicity versus the [Fe/H] as derived 
with the use of Monte Carlo simulations. The circles indicate the 
mean metallicity in each metallicity bin, while the eiTor bars in- 
dicate the standard deviation of the mean. 



are actually old AGB stars ( Durrel l, Harris & Pritchet 120011 
Martinez-Delgado & Apai'icio |1997[ and references therein). In 
order to quantify the effect of the presence of such stars, we con- 
struct the MDF of KDG61 using Padova isochrones (Girardi 
et al. 120081 Marigo et al. 2008 ), which also include the AGB 
phase. We run the interpolation code once using isochrones that 
only include the RGB phase and once with isochrones that only 
include the AGB phase, for a constant age of 12.5 Gyr and a 
range in metallicities from [Fe/H] - -2.36 dex (or Z = 0.0001) 
to [Fe/H] = -0.54 dex (or Z = 0.006), with a step of 0.1 dex 
in [Fe/H]. The derived mean values differ only by 0.04 dex in 
[Fe/H] with a mean of ([Fe/H]) = -1.24 dex for the MDF con- 
structed using isochrones that include only the RGB phase. The 
MDF of the stars that were fit using isochrones that include only 
the AGB phase becomes more metal-rich. Furthermore, the de- 
rived 1 cr spreads in [Fe/H] have comparable values of 0.26 dex 
when we include only the RGB phase and of 0.27 dex when we 
include only the AGB phase. In addition, if we randomly assign 
22 % of the stars within the RGB box with metallicities as de- 
rived using only the AGB phase, while the remaining 78 % of 
the stars with metallicities as derived using only the RGB phase, 
then the resulting MDF has a mean of ([Fe/H]) = -1.24 dex 
with 1 cr spread in [Fe/H] of 0.29 dex. This mean metallicity is 
comparable to the one we compute when we use only the RGB 
phase to derive the metallicities. The shape of the MDFs in all 
these cases does not change significantly. Thus, we can safely 
conclude that the presence of these contaminating old AGB stars 
within the RGB box does not aff'ect the derived MDFs' proper- 
ties significantly. 

In Fig. 12] we overplot the metallicity distribution convolved 
with the eiTors in metallicity (solid line). Also shown in Fig. |2] 
(dashed lines) are fits of Gaussian distributions with the ob- 
served mean and dispersion. For each dSph we compute the 
mean metallicity, ([Fe/H]), as well as the error- weighted mean 
metallicity, ([Fe/H]),t,, along with the corresponding intrinsic 1 cr 
dispersions. We show them in Table |3] columns (2) and (3), 
while the error bars in Fig. |2] indie ate the 1 cr dispersion of the 
error-weighted mean metallicities. The errors in metallicity are 
computed from a set of Monte Carlo simulations, in which each 
star is varied by its photometric uncertainties (both in color and 
magnitude, as given by the Dolphot output) and re-fit using the 



Table 3. Derived Properties. 



Galaxy 


([Fe/H]) ± a 


([Fe/H])H.±tr 


K-S" 


Iagb 




(dex) 


(dex) 


(%)" 




(1) 


(2) 


(3) 


(4) 


(5) 


KDG61 


-1.65 ±0.28 


-1.49 ±0.26 


16 


0.07 


KDG64 


-1.72 ±0.30 


-1.57 ± 0.23 


12 


0.09 


DD071 


-1.64 ±0.29 


-1.56 ±0.24 





0.09 


F12D1 


-1.56 ±0.27 


-1.43 ±0.34 


8 


0.07 


DD0 78 


-1.51 ±0.35 


-1.36 ±0.20 





0.09 


DD044 


-1.77 ±0.29 


-1.67 ±0.19 





0.11 


IKN 


-1.38 ±0.37 


-1.08 ±0.16 




0.08 


F6D1 


-1.63 ±0.30 


-1.48 ±0.43 


0.036 


0.03 


HS 117 


-1.65 ±0.32 


-1.41 ±0.20 


55 


0.14 



" The probabilities indicate whether the populations under conside- 
ration are from the same distribution. 



Table 4. Error-weighted mean metallicities for the 10.5 Gyr and 
8.5 Gyr isochrones MDFs. 



Galaxy 


([Fe/H]),,,,o.5±tr 


([Fe/H]),,8.5 ± 0- 




(dex) 


(dex) 


(1) 


(2) 


(3) 


KDG61 


-1.37 ±0.27 


-1.32 ±0.27 


KDG64 


-1.45 ±0.25 


-1.39 ±0.26 


DD071 


-1.41 ±0.25 


-1.34 ±0.27 


F12D1 


-1.31 ±0.35 


-1.25 ±0.36 


DD0 78 


-1.28 ±0.21 


-1.24 ±0.22 


DD044 


-1.60 ±0.20 


-1.53 ±0.20 


IKN 


-1.03 ±0.17 


-0.98 ±0.17 


F6D1 


-1.37 ±0.48 


-1.31 ±0.50 


HS 117 


-1.31 ±0.21 


-1.27 ±0.22 



identical isochrone interpolation as described above. The 1 cr 
scatter of the output random realizations was then adopted as 
the metallicity eiTor for each star. In the right panel of Fig. |4] 
we show the errors in metallicity computed as described above 
versus the metallicities derived for all the stars within the RGB 
box, here for KDG 64 as an example. The eiTor in metallicity 
increases towards the metal-poor part, which is due to the spac- 
ing between the isochrones that becomes narrower towards the 
metal-poor part. 

In order to further quantify the effect of the assumption of 
the constant age on the MDFs, we apply again the same anal- 
ysis using two different constant ages for the isochrones in the 
interpolation method. The first constant age for the isochrones is 
10.5 Gyr, while the second constant age is 8.5 Gyr. We repeat the 
isochrone interpolation with the bounding boxes and the metal- 
licity ranges being kept the same in all cases. We show the results 
for the MDFs of all the dSphs in Fig.|5] where the MDFs for the 
10.5 Gyr isochrones are shown in the left panels and the ones for 
the 8.5 Gyr isochrones in the middle panels. The derived error- 
weighted mean metallicities ([Fe/H]),vjo.5 and ([Fe/H])vv_8.5, for 
the 10.5 Gyr and 8.5 Gyr isochrones, respectively, are shown in 
TablelH along with their corresponding dispersions. 

In addition, the star-by-star difference in [Fe/H] as derived 
using the 10.5 Gyr and 8.5 Gyr isochrones is shown in Fig.|5] 
right panels. The maximum difference in the derived [Fe/H] us- 
ing the 10.5 Gyr isochrones minus the isochrones with a con- 
stant age of 12.5 Gyr is less than 0.20 dex in all the cases, while 
the maximum difference in the derived [Fe/H] using the 8.5 Gyr 
isochrones minus the isochrones with a constant age of 12.5 Gyr 
is less than 0.40 dex in all the cases. Finally, the overall shape of 
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the MDFs as derived for the 10.5 Gyr and the 8.5 Gyr isochrones 
does not change significantly. 

3.3. Population gradients 

In order to examine the presence or absence of population gradi- 
ents in our dSph sample, we construct the cumulative histograms 
of the stars in each dSph selected in two metallicity ranges, de- 
fined as above and below the respective ([Fe/H])„,. Since the 
dSphs can be considered as being elliptical in projection to first 
order, we define in the following the elliptical radius r as 

where x and y are the distance along the major and minor axis, 
and € is the ellipticity. The major and minor axes are computed 
by fitting an ellipse to contours of the number counts of all stars 
above the 1 cr level. This ellipse is shown in white in Fig. |6] and 
represents the elliptical shape that was chosen for each dSph. In 
the same Fig. |6] we show the contours above the 0.5 cr to 2 cr 
level, which are overlaid on top of density maps. In the study of 
population gradients we exclude the IKN dSph since the field of 
view does not cover the whole extent of the galaxy and further- 
more is contaminated by a bright foreground star. In addition, 
we do not show the elliptical shape for IKN. 

We show the cumulative metallicity distributions in Fig. |7] 
and Fig. |8] (middle panels). We show in the same figures the ra- 
dial metallicity distributions (upper panels) and the mean radial 
metallicity profiles (lower panels). Each radial metallicity profile 
shows the mean values of metallicity within an elliptical annu- 
lus versus the elliptical annulus in units of the effective radius 
r^ff. The values for r^.// are listed in the column (10) of Table|2] 
The error bars in the metallicity profile correspond to the stan- 
dard deviation of the mean metallicity in each elliptical radius 
annulus. 

3.4. Density Maps 

We now examine the spatial distribution of the stellar popula- 
tions, separated into a metal-poor and metal-rich component. 
For that purpose, we define two stellar populations, the first in- 
cludes stars having a metallicity less than or equal to the value of 
-1.80 dex ("metal-poor"), while the second includes stars with 
a metallicity larger than or equal to -1.30 dex ("metal-rich"). 
All galaxies have peak values that lie well in between those cuts 
so that the metal-poor and the metal-rich tails are representa- 
tively sampled for all dSphs. For these two populations we con- 
struct the gaussian-smoothed density maps, shown in Fig.|9]and 
Fig-Hni upper and middle panels for each dSph. 

3.5. Luminous AGB stars 

As already noted before, all the dSphs in our sample contain a 
significant number of luminous AGB star candidates. These stars 
are located above the TRGB and have ages ranging from 1 Gyr 
up to less than 10 Gyr. We broadly refer to stellar populations in 
this age range as "intermediate-age" populations. Assuming that 
the metallicities of dwarf galaxies increase with time as star for- 
mation continues, we may also assume that these intermediate- 
age populations are more metal-rich than the old populations. 
We note, however, that dwarf galaxies do not necessarily expe- 
rience smooth metal enrichment as a function of time (see, e.g., 
Koch et al. ,2()07atl2007bl l. 



In Fig. |9] and Fig. [10] the density maps in the lower pan- 
els show the spatial distribution of these intermediate-age stars 
for each dSph. We consider as luminous AGB stars the stars 
that are brighter by 0.15 mag than the TRGB (ArmandrofF et 
al. 1 19931 1 and that lie within 1 mag above {Jtrgb - 0.15) mag. 
In addition, we consider stars within the color range of a < 
(y - I)q < a + 2.50 (mag), where the left-hand limit a is equal 
to the color of the TRGB of the most metal-poor isochrone we 
use, dereddened for each dSph using the extinction values listed 
in columns (6) and (7) of Table |2] Then, the right-hand limit is 
the left-hand limit plus 2.50 mag. This selection criterion was 
motivated by the work of Brewer, Richer & Crabtree (|1995i l and 
Reid & Mould (fT984l i. 



4. Discussion 

4.1. Photometric Metaliicity Distribution Functions 

The photometric MDFs in Fig. |2] indicate that these dSphs cover 
a wide range in metallicity. All of them seem to have a steeper 
cut-off in their metal-rich end. This can be easily seen if we com- 
pare the MDFs to the fitted gaussian distributions, indicated by 
the dashed lines in Fig.|2l We do not expect the MDFs to follow 
a gaussian distribution since they are shaped by the star forma- 
tion histories of each dSph. For instance, a steep cut-off of the 
MDF toward the metal-rich tail could be indicative of the oc- 
currence of strong and continuous galactic winds (Lanfranchi 
& Matteucci 2003j [20071 Koch et al. [2006l l or of the effects of 
extended star formation and SNe la "pockets" of localized, in- 
homogeneous enrichment (Marcolini et al. ' 20081 ). 

The low mean metallicities, ([Fe/H]), that are derived from 
the distribution functions and are shown in Table [3] columns 
(2) and (3), indicate that the M81 dSphs are metal-poor sys- 
tems, which points to a low star formation efficiency in analogy 
with the LG dSphs (e .g., La nfran chi & M atteucci 2003. Grebel, 
Gallagher & Harbeck [2U03l Koch [2009l l. 

One exception is the dSph IKN, which shows a high mean 
metallicity for its luminosity. Objects that have high metallic- 
ity for their luminosity and that, most importantly, are dark 
matter free are promising candidates for tidal dwarf galaxies. 
Their properties are set by their formation mechanism. They 
are believed to form out of the dark matter free material that 
was expelled during the tidal interaction of the parent galax- 
ies (Bournaud et al. 120071 and references therein). One possi- 
ble tidal dwarf galaxy candidate has been identified in the M 8 1 
group, namely HolmberglX (Makarova et al. 120021 Sabbi et 
al. 2008) and these are favoured to be detected in such recently 
interacting groups. These systems contain a young stellar com- 
ponent, while their older stellar populations are believed to con- 
sist of stars from their parent galaxies, which is M81 in the 
case of Holmberg IX. There is no information available about 
the presence or absence of dark matter in this system. 

In the case of IKN, we should consider the fact that its 
stellar metallicity bears the imprint of the medium that formed 
these old stars, while young stars are not observed in this dwarf. 
That makes it distinct from young tidal dwarf candidates like 
Holmberg IX. A connection with the recent interactions of the 
M81 group is not obvious. IKN might be an old tidal dwarf 
galaxy if such systems exist. Alternatively, IKN may have un- 
dergone substantial mass loss in the past, leaving it as a low- 
luminosity but comparatively high-metallicity dSph. Without 
data on its detailed structure and kinematics, we cannot distin- 
guish between these possibilities. 
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Fig. 6. Contour plots, shown in red, for the nine dSphs which are overlaid on top of density maps. The elliptical shape chosen 
for each dSph, with the exception of IKN, is shown with the white ellipse. The star symbols in some plots coiTespond to bright 
foreground stars, while in the case of KDG64 corresponds to a background galaxy. The diamond symbol in the case of F6D1 
corresponds to an extended background galaxy. The unit of the colorbars is number of stars per (50 pixels)^. 



The metallicity spreads of the studied dSphs are large, span- 
ning 1 o" ranges from 0.27 dex to 0.37 dex, or intrinsic, eiTor- 
weighted 1 cr ranges from 0.16 dex to 0.43 dex. These abun- 
dance spreads are comparable to the ones observed in the LG 
dSphs (Grebel, Gallagher & Harbeck l2003l Koch 1200 9 ) and 
may indicate the presence of multiple stellar populations and / or 
extended star formation histories. According to the models of 
Mai-coHni et al. (2008 ) and Ikuta & Arimoto (2002 ), the initial 
star formation in dSphs may have lasted as long as 3 Gyr or 
even longer, which can lead to a large dispersion in [Fe/H]. For 
ages older than 10 Gyr, the shape of the MDF does not depend 
strongly on age as described in Sec. l3.2l and shown in Fig.|5] 

4.2. Luminosity-Metallicity Relation 

In Fig. [TT] upper panel, we show the luminosity-metallicity re- 
lation compiled for the dwarf galaxies in our LG as studied by 
Grebel, Gallagher & Harbeck (120031 ). with the addition of thir- 
teen dSphs in the M 8 1 group. This compilation includes LG ob- 



jects with mean metallicities derived from either spectroscopic 
or photometric studies. The mean metallicities for the nine M 8 1 
group dSphs listed in Table [3j column (3), are from this work, 
while the metallicities for the remaining four dSphs are adopted 
from the literature and are computed using the mean (V - I)o 
color of the RGB stars at the luminosity of M/ = -3.5 mag (from 
Caldwell et al. ( 119981 ) for BK5N and F8D1 and from Sharina et 
al. (2008 ) for KKH 57 and BK6N). 

Overall, the M81 group dwarfs follow the luminosity- 
metallicity relation quite well, albeit some of them exhibit a ten- 
dency of being slightly more metal -poor than LG dSphs of com- 
parable luminosity. Therefore, they mainly populate the region 
defined by the LG dSphs while a few are located in the border re- 
gion between the dSph and dirr locus defined by the LG dwarfs. 
The M 8 1 group dSphs that seem to lie in this apparent transition 
region are KDG 6 1 , KDG 64, DDO 44 and DDO 7 1 . Out of these 
four objects, three are classified as transition-types, namely, 
KDG 61, K DG 64 and DDO 71 (Karachentsev & Kaisin 120071 
Boyce et al. l2001l ) based on HI detections and Ha emission. Also 
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Fig. 7. In each panel from left to right and from top to bottom we show the radial metallicity distributions (top panels), their 
cumulative distributions (middle panels), and the radial mean metallicity profile (bottom panels). The cumulative histogram is for 
stars selected in metallicity above and below the weighted mean value ([Fe/H])^^. listed in Table [3] 
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Fig. 8. Same as in Fig.|2]for the remaining two dSphs. Note that IKN is excluded from this analysis. 



among the dwarfs that coincide with the LG dSph locus, one 
dwarf is classified as transition-type, namely HS 117, with HI 
associated with it (Huchtmeier & Skillman 119981 Karachentsev 
et al. 12006) . 

Transition-type dwarfs are galaxies that share properties of 
both morphological types. Their stellar populations and star for- 
mation histories resemble those of dSphs and their gas con- 
tent and present-day star formation activity is akin to low-mass 
dlrrs. It has been suggested that transition-type dwarfs are indeed 
evolving from dlrrs to gas-deficient dSphs. 

The projected spatial distribution of the dSphs within the 
M81 group is shown in Fig. 1 of Karachentsev et al. (I2002I I. 
while their three-dimensional view is shown in their Fig. 6. In 
Fig. [TT] middle panel, we plot the deprojected distances from 
the M 8 1 galaxy, R, versus the mean metallicities for the 1 3 M 8 1 
group dSphs. The deprojected distances from the M81 galaxy, 
R, are adopted from Karachentsev et al. (I2002I I and are listed in 
Table|2] column (9). The most distant dSph is DDO 44 which be- 
longs to the NGC 2403 subgroup, while KDG 61 is the one clos- 
est to the M81 itself. Interestingely, the dwarf KDG 61 which 
is classified as a transition-type based on HI detections, is the 
closest to the M 81 itself, with a deprojected distance of 44 kpc 
(Karachentsev et al. I2002I I. The remaining three dwarfs classi- 
fied as transition-types, namely KDG 64, DDO 71 and HS 117, 
lie in a deprojected distance of more than 100 kpc. As discussed 
already, the most metal-rich dSph studied by us is IKN, while 
according to the value that Caldwell et al. (1998 ) provide for 
F8D1, this is the most metal-rich dSph in this group so far stud- 
ied. We see no trend of the mean metallicity with the deprojected 
distance. 



4.3. Population Gradients 

By examining the cumulative metallicity distributions in the 
middle panels in Fig. [T] and Fig. [8] we conclude that the metal- 
licity gradients are present in the case of DDO 71, DDO 78, 
DDO 44 and F6D1, while the metallicity gradients are less 
pronounced or not present in the remaining dSphs. We again 
separate the RGB stars in each dSph into two samples, 
where we choose to separate the distributions at the observed 
weighted mean metallicity. The probabilities from the two-sided 
Kolmogorov-Smirnov (K-S) test that the two components are 
drawn from the same parent distribution are listed in Tabled col- 
umn (4). The K-S results are consistent with showing spatially 
separated populations in the case of DDO 71, DDO 78, DDO 44 
and F6D1. In the case of the remaining dSphs, the gradients are 
less pronounced and the metal-rich and metal-poor populations 
have different distributions at the 84 - 99.7 % confidence level 
(>1.5 cr), except for HS 117. In all cases, in which we observe 
such metallicity segregation, the sense is that more metal rich 
stars are more centrally concentrated, as also found in the ma- 
jority of th e LG d Sphs (Harbeck et al. 120011 Tolstoy et al. l2004l 
Ibata etal. l2006l l. 

4.4. Density Maps 

The density maps in Fig. |9]and Fig. [10] are useful to study the 
spatial distribution of the metal-rich (upper panels) and metal- 
poor (middle panels) population of each dSph. From these den- 
sity maps we conclude that each dSph has a different stellar spa- 
tial distribution of their metal-rich and metal-poor stellar com- 
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ponent. All of them show either a spatial variation of the cen- 
troids of the two stellar populations, as is the case of F12D1 
and KDG 64, or that the metal-rich population is more centrally 
concentrated, as is the case of DD071, DD044 and F6D1. 
These findings agree well with the ones from the cumulative his- 
tograms, though we should keep in mind that the metal-poor and 
metal-rich stellar populations involved are differently selected. 
DD078 and IKN are clearly fairly metal-rich while KDG 61, 
KDG 64 and DDO 44 have prominent metal-poor populations. 

4.5. Luminous AGB stars 

Luminous AGB stars were also detected in two more dSphs in 
the M81 group, namely BK5N and F8D1 (Caldwell et al. 1998). 
These luminous AGB stars may include carbon stars and may 
be long-period variables (LPVs) as have been found in other 
early-type dwarf galaxies (e.g., Menzies et al. 120021 Rejkuba et 
al. 120061 Whitelock et al. '2009'), but we can not establish this 
for certain based on our current data. We compute the fraction 
of the luminous AGB stars, /acb, defined as the number of the 
luminous AGB stars, Nacb, counted within the magnitude bin 
considered in Sec. 13.51 over the number of the RGB stars within 
one magnitude below the TRGB, Nrcb- In order to estimate the 
Nrcb, we take into account that approximately 22 % of the stars 
we count within one magnitude below the TRGB are old AGB 
stars (DuiTell, HaiTis & Pritchet 2001). Thus, the Nrcb is equal 
to 78 % of the stars we count within one mag below the TRGB. 
The fractions /acb we derive in this way are listed in Table [3] 
column (5). 

4.5.1 . Blends, blue straggler progeny and foreground 
contamination 

We now discuss the possibility that these luminous AGB stars 
may actually be (1) blends of bright RGB stars (Renzini ll998l) . 
(2) blue straggler progeny (Guarnieri, Renzini & Ortolani[1997J 
and references therein), or (3) due to the foreground contamina- 
tion. 

In order to evaluate case (1), we use our artificial star exper- 
iments in order to quantify the number of the blends of bright 
RGB stars, Nbknds, that may contribute to the detected number 
of the observed luminous AGB stars. We only consider here the 
case of a blend of two equal-magnitude stars. The magnitude 
of the blended star is always 0.75 mag brighter than the initial 
magnitude of the two superimposed stars. 

We want to determine the location in the CMD of all the 
RGB stars that can end up as blends within the location in the 
CMD of the luminous AGB stars, as defined in Sec. l3.5l and fur- 
ther called as luminous AGB box. For that purpose, if we assume 
that the stars within the luminous AGB box were all blends, then 
they would originate from stars that have magnitudes 0.75 mag 
fainter. Thus, we shift the luminous AGB box by 0.75 mag to- 
wards the fainter magnitudes and furthermore we only consider 
the stars with magnitudes fainter than the TRGB. This procedure 
defines the location of the RGB stars that can end up as blends 
within the luminous AGB box and we call this the "RGB blends 
box". 

From the stellar catalogue we use as an input for the arti- 
ficial star experiments, we select the stars that have such input 
magnitudes to place them within the "RGB blends box". From 
these input stars, we consider as blends the ones that have output 
magnitudes that can place them above the TRGB. We normal- 
ize the number of these blends to the number of the total input 



stars that are located within the "RGB blends box". Finally, the 
number of the observed blends is proportional to the number of 
the observed RGB stars located within the same "RGB blends 
box". Thus, the Nbiends for all the dSphs computed this way is 
equal to 5 blends, 1 1 blends, 12 blends and 2 blends in the case 
of DDO 44, DDO 78, IKN and HS 117, respectively, while in all 
the other cases the number of blends is less than 1 . Thus, the 
fraction of the blends defined as the number of blends divided 
by the number of the RGB stars within 1 mag below the TRGB, 
is less than 0.6 % in all cases but for IKN which is equal to 0.9 %. 

In the case (2), Guarnieri, Renzini & Ortolani (1997) point 
out that the number of blue straggler progeny is of the order of 
~2 % of all stars that reach the luminous AGB phase. 

In the case (3), we estimate the foreground contamination 
using the TRILEGAL code (Vanhollebeke, Groenewegen & 
Girardi l2009t Girardi et al. 120051 ). We count the number of fore- 
ground stars that fall within the same location in the CMD as in 
the luminous AGB box and these are considered to be the num- 
ber of expected foreground contamination. In all the cases, the 
number of foreground stars is 4, with the exception of DDO 7 1 
and DDO 44 where the number of foreground stars is 3 and 5, 
respectively. This translates to a fraction of foreground stars, de- 
fined as the number of foreground stars divided by the number 
of RGB stars within 1 mag below the TRGB, of less than 0.7 %, 
with the exception of F6D1 where this fraction is equal to 2 %. 

We do not consider the case of old AGB LPVs, whose large 
amplitude variations may place them above the TRGB (Caldwell 
et al. 11998b . as an additional source of contamination in the lu- 
minous AGB box, since the studied dSphs have mean metallici- 
ties of less than -1 dex. Such old AGB LPVs, with ages greater 
than or equal to 10 Gyr, have been observed above the TRGB in 
metal-rich ([Fe/H] > -1 dex) globular clusters (e.g., Guarnieri, 
Renzini & Ortolani [19971 and references therein). 

We can now add all the contributions estimated in the above 
three cases, for each dSph. We call the sum of these three contri- 
butions the number of total contaminants, N com. mi, and compute 
their fraction f com, tot - N com, tot I Nrcb- The fraction of the to- 
tal contaminants is less than 1 % in all cases apart from IKN, 
HS 117 and F6D1, where the fraction of the total contaminants 
is approximately 1.1 %, 1.3 % and 2.1 %, respectively. We note 
that in all cases, there is a significant fraction of luminous AGB 
stars that can not be accounted for by considering the contribu- 
tion of blends, binaries and foreground contamination. The dSph 
F6D1 is an exception to that, where the fcom.tot is ~ 2 %, as com- 
pared to the /agb which is ~ 3 %. We note though that in the case 
of F6D 1 the number of stars counted in the luminous AGB box is 
equal to 6 stars. Thus, we conclude that the luminous AGB stars 
are a genuine population for the majority of the dSphs studied 
here. 

4.5.2. Luminous AGB density maps and fractions 

From the density maps of the luminous AGB stars shown in 
Fig. |9] and Fig. [10] lower panels, we see that if we consider the 
peak densities or the bulk of the luminous AGB stars, then it 
seems that these are more confined to the central regions of the 
dwarfs, a behaviour similar to what is found for the Fornax dSph 
(Stetson, Hesser & Smecker-Hane 119981 ). If we consider the 
overall distribution then we note that for most of the dSphs these 
stars are rather more widely distributed, following the distribu- 
tion of the metal-poor stars, with the exception of KDG 64 and 
DDO 71 where their AGB stars' distributions coincide mostly 
with the metal-rich population, which in the case of DDO 7 1 is 
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centrally concentrated. We conclude that the intermediate-age 
stellar component is well-mixed with the old stellar component. 

This behaviour is similar to the AGB stars' spatial distribu- 
tion of the LG dwarfs. Indeed, Battinelli & Demers (I2004at and 
references therein) discuss that for the LG dwarfs, for which 
there are carbon star studies, their Carbon-rich stars are dis- 
tributed such that they coincide with the spatial distribution of 
the old stellar component. An exception is the dE NGC185, 
where the AGB stars are concentrated more in the centre than 
the old stellar component (Battinelli & Demers l2004bl l. a simi- 
lar behaviour as observed in the two M 8 1 group dSphs discussed 
above. 

We plot the fACB versus the deprojected distance from M 81, 
R, in the lower panel of Fig[TT] The highest fraction of luminous 
AGB stars is observed in HS 117 and the lowest one in KDG61 
and F12D1. We do not see any trend of the /agb with increasing 
deprojected distance from M 8 1 . If we compute the net fraction 
of the luminous AGB stars, by subtracting the fraction fcont,tot, 
due to the contribution of blends, binaries and foreground con- 
tamination, from the fraction fACB hsted in the column (5) of 
Table [3] none of the trends and conclusions change. 



5. Conclusions 

We use the CMDs of nine dSphs in the M 8 1 group to construct 
their photometric MDFs. These MDFs show populations cov- 
ering a wide range in metallicity with low mean metallicities 
indicating that these are metal-poor systems. All MDFs show 
a steeper fall-off at their high-metallicity end than toward their 
low-metallicity end indicating that galactic winds may play a 
role in shaping their distribution. 

We compute the mean metallicity, ([Fe/H]), and the mean 
metallicity weighted by the metallicity error, ([Fe/H])v^., along 
with their coiTesponding standard deviations. The most metal- 
rich dSph in our sample is IKN even though it is the least lumi- 
nous galaxy in our sample. IKN's comparatively high metallic- 
ity may indicate that it is a tidal dwarf galaxy or that it suffered 
substantial mass loss in the past. We do not see any correlation 
between the ([Fe/H]) and the deprojected distance from the M 8 1 
galaxy, R. 

We use the mean metallicity weighted by the metallicity er- 
rors, ([Fe/H])n , to select two stellar populations having metal- 
licities above and below that value. For these two stellar popu- 
lations we construct cumulative histograms, as a way to search 
for population gradients in metallicity. We find that some dSphs 
show strong metallicity gradients, while others do not. In dSphs 
with radial metallicity gradients the more metal-rich populations 
are more centrally concentrated. 

Furthermore, we study the spatial (i.e., two-dimensional) 
distribution of our defined metal-rich and metal-poor stellar pop- 
ulations. This refined look no longer assumes radial symmetry, 
and we now find that in some dwarfs the metal-rich population 
is more centrally concentrated, while others show offsets in the 
centroid of the two populations. By examining the distribution 
of the luminous AGB stars, we conclude that, for the majority of 
the dSphs, these stars have mostly extended distributions, indi- 
cating that they have been well-mixed with the metal-poor stellar 
population. We do not find any correlation between the fraction 
of luminous AGB stars and the deprojected distance from the 
M 8 1 galaxy. While present-day distances may not be indicative 
of the dwarfs' position in the past and while their orbits are un- 
known, the apparent lack of a correlation between distance and 
evolutionary history may suggest that the evolution of the dwarfs 



was determined to a large extend by their internal properties and 
not so much by their environment. 

Finally, there are some M81 dSphs that straddle the 
transition region between LG dSphs and dirrs in the 
metallicity-luminosity relation. We may be observing low- 
luminosity transition-type dwarfs moving toward the dSph lo- 
cus. Interestingly, these dwarfs are slightly more luminous than 
the bulk of the LG transition dwarfs. Perhaps some of the M 81 
dwarfs experienced gas stripping during the recent interactions 
between the dominant galaxies in the M 8 1 group. 
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Fig. 5. Metallicity distribution functions for the nine dSphs com- 
puted in the same way using isochrones with two different ages: 
left panels for a constant age of 10.5 Gyr and middle panels for 
a constant age of 8.5 Gyr. The solid lines show the metallicity 
distribution convolved with the errors in metallicity. The error 
bars correspond to the 1 cr spread for the weighted mean metal- 
licity we derive from our data. Right panels: Star-by-star differ- 
ence of the derived metallicities, A[Fe/H], using the 10.5 Gyr 
isochrones minus the 12.5 Gyr isochrones versus the [Fe/H] of 
the 10.5 Gyr isochrones, indicated with the open circles. The 
star-by-star differences of the 8.5 Gyr isochrones minus the 
12.5 Gyr isochrones are indicated with the dots. 
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Fig. 11. Upper. Luminosity-metallicity relation for LG dwarf 
galaxies, after Grebel, Gallagher & Harbeck (2003 ), together 
with the 13 dSphs of the M81 group. LG dSphs are plotted 
with blue asterisks, LG dirrs are shown as green crosses, and red 
dots indicate the available M 8 1 data. Nine out of the thirteen 
M81 group dSphs have been studied here, while the remain- 
ing four, marked with a red circled dot, have been adopted from 
the literature, as discussed in the text. The red squared dots and 
green squared crosses indicate the transition-types of the M 8 1 
group and LG, respectively. Middle: Mean metallicities versus 
the deprojected distance from M81, R, for the 13 M81 group 
dSphs. The circled dots correspond to the four dSphs for which 
the metallicities were adopted from the literature, as discussed 
in the text. Lower. Fraction of AGB stars versus the RGB stars 
within 1 mag below the TRGB, fACB, versus the deprojected dis- 
tance from M 8 1 , R, for the nine M 8 1 group dSphs studied here. 
With the circled dot we show the corresponding fraction /acb in 
the case of F6D 1 . 
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Fig. 9. From left to right and from top to bottom we show the density maps for each dSph, after smoothing with a Gaussian kernel. 
In each upper and middle panel, the "metal-rich" population corresponds to stars having [Fe/H]> -1.30 dex while "metal-poor" 
refers to [Fe/H]< -1.80. The bottom panel corersponds to the density map of the luminous AGB stars, as defined in the text. The 
unit of the colorbars is number of stars per (100 pixels)^. 
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Fig. 10. The same as in Fig.|9]but for the remaining three dSphs. Note that the zero density region in the center of IKN is an artifact 
due to a bright star contaminating the field of view. 
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